SISSA 60/2004/FM 



in 
o 
o 

(N 
D 

(N 
(N 

(N 
> 

o 



^ : 
o . 

Oh: 

> 



Form factors in the massless coset models 

su{2)k+i su{2)k/su{2)2k+i 
Part I 

Paolo Grinza"'^ ^ and Benedicte Ponsot"'"^ ^ 

""International School for Advanced Studies (SISSA), 
Via Beirut 2-4, 34014 Trieste, Italy 
INFN sezione di Trieste 

Laboratoire de Physique Mathematique, Universite Montpellier II, 
Place Eugene Bataillon, 34095 Montpellier Cedex 05, France 

'^Department of Theoretical Physics, 
Research School of Physical Sciences and Engineering, 
Australian National University, 
Canberra, ACT 0200, Australia 



Massless flows between the coset model su{2)k+i ® su{2)k / su{2)2k+i and the minimal 
model Mk+2 are studied from the viewpoint of form factors. These flows include in particular 
the flow between the Tricritical Ising model and the Ising model. Form factors of the trace 
operator with an arbitrary number of particles are constructed, and numerical checks on 
the central charge are performed with four particles contribution. Large discrepancies with 
respect to the exact results are observed in most cases. 




Abstract 



PACS: ll.lO.-z ll.lO.Kk 11.55.Ds 



grinza@lpm.univ-montp2.fr 
^ ponsot@fm.sissa.it, benedicte.ponsot@anu.edu.au 



1 



Introduction 



The general problem of finding a suitable description of Renormalization Group (RG) flows 
between different non-trivial critical point in Quantum Field Theory is an open and appealing 
problem. In this respect 2d integrable QFTs provide a privileged framework where there is the 
actual possibility of finding a detailed description of such flows. As a matter of fact, one of the 
remarkable consequences of integrability is the knowledge of the exact S'-matrix which allows 
to calculate exact form-factors [1-3]. From the latter one can reconstruct correlators by means 
of the spectral expansion which provides a non-perturbative representation of them. Such a 
representation has often been very accurate in the past -though it was noticed some time ago 
in [4] that this common belief, even for massive theories, is too optimistic in general- and allows, 
for example, to recover the conformal data using the c-theorem sum rule [5, 6] which gives the 
difference between the central charges of the UV and IR fixed points. Integrable RG flows which 
limit in the infrared is a non-trivial CFT can be described by massless excitations, for which an 
exact scattering matrix can be found. More precisely, we deal with right and left movers and 
three different types of scatterings, associated to right-right, left-left, and right-left interactions. 
The right-right and left-left S'-matrices, whose definition is formal, are independent of the RG 
scale, and are solely characterized by the properties of the IR fixed point CFT. On the contrary, 
the right-left scattering is quite rigorously defined [7]. It becomes trivial in the IR limit, thus the 
left and right movers decouple, and we obtain the IR CFT. In [8], using the S'-matrix proposed 
in [7] , a few form- factors^ of some local and non local operators were constructed for the simplest 
model describing the flow [9-11] between the Tricritical Ising model (TIM) and the Ising model 
(IM), where the supersymmetry is spontaneously broken. Remarkably enough, the numerical 
results performed in [8] show that the knowledge of the form-factor with the lowest number 
of intermediate particles is enough to give quite an accurate approximation for the correlation 
function at almost all RG scale. Of course it is desirable to figure out whether these spectacular 
features remain valid in less trivial massless flows. In this respect, one has to take into account 
results obtained recently in [12] about the calculation of form factors in a massive integrable 
model of QFT called the SS model [13], as well as in its RSOS restrictions. The numerical 
checks performed in [12] show in particular the two-particles approximation of the c-theorem 
sum rule fails to give the usual good approximation of the exact results, showing instead large 
discrepancies (about 20-25%, see below). 

In this article we will consider the construction of form factors of the trace operator for an 
arbitrary number of intermediate particles in the family of massless flows [14] between the UV 
coset models [15] su(2)fc+i (g) su{2)f:/su{2)2k+i and the minimal models Mfc_|_2 in the IR. These 
flows, whose direction is given by the irrelevant operator TT, include in particular for k = 1 the 
above mentioned flow from TIM to IM, and for k = -\-oo the massless flow from the Principal 
Chiral Model at level 1 (PCMi) to the SU {2)i WZNW model [16]. 

We recall that the form factors of the trace operator were actually first obtained in [17] in a 
different representation from the one presented in this article; the representation proposed here 
is simpler. Let us mention also that no numerical checks were performed in [17]. 

This article is organized as follows: in section ^ we recall known facts on the construction of 
form factors in the Sine-Gordon model and its RSOS restriction, that we will need in section |21 
where we construct the form factors of the trace operator in the massless flows mentioned above. 
Finally, in sect. |21 we perform checks on the four particle form factor by making a numerical 

■^We refer the reader to [8] for a discussion on form-factors in massless theories. 
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estimation of the variation of the central charge between the UV and the IR CFTs; we will see 
that the accuracy becomes poorer and poorer as one increases k, in an even worst manner than 
what was first observed in the massive case [12], and we will provide a reasonable explanation 
for such a loss of accuracy. 



1 Form factors in the sine-Gordon model and RSOS restriction 

In this section we recapitulate known results on form factors in the SG model in the repulsive 

regime. 

The Sine-Gordon model alias the massive Thirring model is defined by the Lagrangians: 



C^™ = i,{i^d-M)i,^-g{i^^,i,)\ (1) 



1 

T 

respectively. The fermi field -0 correspond to the soliton and antisoliton and the bosc field ip to 
the lowest breather which is the lowest soliton antisoliton bound state. The relation between 
the coupling constants was found in [18] within the framework of perturbation theory: 

_ _ TT 

The two particles S-matrix of the SG model [19] contains the following scattering amplitudes: 
the two-soliton amplitude ap{0), the forward and backward soliton anti-soliton amplitudes bp{9) 
and Cp{9): 

, fn\ sinhe/p sinhivr/p 

smh(i7r — u)/p smh(i7r — u)ip 



ap{e) = exp / 

^0 



°° disinhi(l -p)tsinhf^ 



t cosh I sinh ^pt 

This S'-matrix satisfies the Yang-Baxter equation as well as the unitarity condition: 

s^p''{e)s^p''{-e) = 1, 

which can be rewritten for the amplitudes as: 

ap{e)ap{-e) = 1, bpie)bp{-e) + cp{e)cp{-e) = i. 

The crossing symmetry condition reads for the amplitudes: 

ap{iTT - 6) = bp{e), Cp{i7r - 0) = Cp{0). 

The repulsive regime (with no bound states) corresponds to the condition p > 1. 

The form factors f{9i, • • • , 02n) of a local operator^ in the SG model are covector valued functions 

that satisfy a system of equations [3], which consist of a Riemann-Hilbert problem: 

fi^l, • ■ ■ di+i, ■ ■ ■ , 92n)Sp'^{0i — 9i+l) = f{Ol, ■ ■■ , 0i+i,9i, ■■■ , 92n), 
f{9l, ■ ■ ■ , 92n-l,92n + 2i7r) = f{92n, 9l, - ■ ■ , 92n-l), 

^Horc wc shall exemplify form factors with an even number of particles only. For form factors with an odd 
number of particles, we refer the reader to [20]. 
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and a residue equation at 9i = ^2n + i^i" '■ 

(2n-l \ 
1 - JJ - &2n) ) eo, eo = Si (g) S2n + Si ® S2n, 

where s (or +) corresponds to the solitonic state (highest weight state), and s (or -) to the 
antisohtonic state. The equations written above do not refer to any particular local operator. 
It is one of the difficulties of the form factor approach to identify solutions of these equations 
with specific operators. 

Form factors containing an arbitrary number of particles for the energy momentum tensor 
were first constructed in [3] by Smirnov. Below, we make the choice to present the posterior 
construction presented in [20,21] by Babujian and Karowski et. al. We refer the reader to these 
two references for more details. 

We first introduce the minimal form factor fp{9i2) of the SG model: it satisfies the relation 

fp{&) = -fp{-OK{e) = U{2m - 6), 

and reads explicitly 

frm • • ^. ^ • • V. ^ sinhi(l 1 - cosht(l - S^) 
fvi^) = -1 smh -flP' (9) = -1 smh - exp / ^ ^ — — f 

Its asymptotic behaviour when 9 ±00 is given by fp{9) Cp €■ ^ ^ , with. tli6 constcint 

dt ( s\i\\\\{l -p)t l-p\ 



. 1 1 r 



t \ sinh^ptcosh^tsinhi pt I 
It is proposed in [20] that form factors in SG can be written^: 

f{9i, ...,9n) = Nnll fp{9ij) f dui... f du^ hp{9, u)pni9, u)^P{9, u), (2) 

•^Cfl «^Cfl 



ICg JCg 

where we introduced the scalar function (completely determined by the S'-matrix) 

2n m 

hp{9,u) =YYY\^4)p{9i-Uj) JJ Tp(Ur-Us), 
i=lj=l l<r<s<m 



with 



fp{u)fp{u + i7r)' ^ (t>p{u)ct>p{~u)' 
^'P{9,u) is the Bethe ansatz state covector: we first define the monodromy matrix Tp as 



A{9i,...,9n.u) B{ 
Ci9i,...,9n,u) D{ 



h\\'.'.\lZu) ) ^Tpi9u...,9n,u)=S^:^'{9,-u)...S^''{9n-u), 



^This representation holds whether operators are local or not, topologically neutral or not. 
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the definition of the Bethe ansatz covector is given by 

"$^{6, u) = 0,l,,,n C{9i, . . . ,9n, ui) . . . C{9i, . . . , 
where Qi,,,n is the pseudo vacuum consisting only of sohtons 

^l...n = S ■ ■ ■ S. 

The number of integration variables m is related to the topological charge q {q £ 7^) of the 
operator considered and the number n of particles through the relation 

q = n — 2m. (3) 

For example, for n = 2 and q = 0, 

^P{9ue2,u) = ^l_{9u92,u)+^'L+{9i,92,u) = b{9i-u)c{92-u)si^S2+c{Gi-u)a{92-u)si®S2. 

It is important to have in mind that the function pn{9,u) is the only ingredient in formula (jSJ 
which depends on the operator considered. Different operators will differ by their p-function. 
If the operator is chargeless, the form factors contain an even number of particles, and if in 
addition the operator is local, then the function satisfies the conditions®: 

1. P2n{9, u) is a polynomial in e^'^^ , (j = 1, . . . , n) and P2n{9, u) = p2n{- ■ ■ ,9i- 2i7r, ... ,u) 

2. P2n{9l = 6'2n+ivr, . . . 92n]Ui . . . Un = 92n) = P2n-2{92 ■ ■ ■ 92n~l]Ui . . . n„_i (^2 • • • 6'2n-l), 
P2n{9l = 6'2n+ivr, . . . 92n]Ui . . . Un = 6'2n±i7r) = P2n-2{92 ■ • .^2n-i;'Ul • ■ ■ Un-l)+P±{92 • • • 6'2„_l), 

where p^''^{92 ■ • • ^2ra-i) are independent of the integration variables. 
3- P2n{9,u) is symmetric with respect to the 0's and the u's. 

4. P2n{9 + In A, u + In A) = A*^ P2n{9, u) where s is the Lorentz spin of the operator. 

Finally, the integration contours Cg consist of several pieces for all integration variables Uj : a 
line from — oo to oo avoiding all poles such that Im^j — vr — e < Imuj < Im^j — vr, and clockwise 
oriented circles around the poles (of the — Uj)) at 9i = Uj, (j = 1, . . . , m). 

Trace of the energy momentum tensor. The trace operator is a spinless and chargeless local 
operator. Its ]5-function is [21]: 

(2n n 2n n \ 

i=i j=i i=i j=i ) 

The residue equation gives the following relation for the normalization Nfn (see [20]): 

^2n - ^2n-2 ^ ^2n " ^2 ^, 4^ • (5) 



®We consider here only the case where the operator is of bosonic type and the particles are of fermionic type. 
If both are fermionic, there is an extra statistic factor to be taken into account, see [21]. 
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The two particles form factor can be computed explicitly: 

2ypjY0 cosh — 

fsci^u) = -^fpiOu) . , 1 \ A si ® S2 + Sl S2), 

smh 2^(i7r - 6'i2) 

in agreement with the result first obtained by diagonalization of the S'^'^-matrix in [1]. The 
normalization for two particles is chosen to be N2 = ^M'^Cp (M being the mass of the soliton), 
in order to have: 

/|g((9i + ivr,6'i) = l-KM^isi O §2 + si 52)- 

Let us note that the two and four particles form factors were checked again Feynman graph 
expansion in [21]. 

Holomorphic components of the stress- energy tensor T±±. It is a chargeless operator with ±2 
Lorentz spin. Its p-function reads [21]: 



-1 



P2n 



E-^N E-^N PsoiO,n). (6) 



\i=i / \j=i / 
This relation ensures the conservation laws: 

5+e = 9_T++, d-e = d+T-~. 

The two particles form factors read: 

fit (0,2) = -^/(gi2)e^(^-+^-) . "Tr"^^ A si .-2 + .-1 s,), 
S smh 2^(i7r - 6'i2) 

and = -^M^C^ = -iVf. 

RSOS restriction [22, 23]. The RSOS restriction describes the <I>i^3-perturbations of minimal 
models of CFT [24] for rational values of p. In particular, we remind that S^^^^ = — 1 is the 
Ising 5-matrix. Form factors in this model can be directly obtained from those of the SG model, 
as explained in [23]. For the trace operator, the RSOS procedure consists of 'taking the half of 
the p- function^ Q), such that its p-function reads: 

2n n 
i=l m=l 

then we should modify the Bethe ansatz state: 

2n 

Ke....e,AO,u) ^ e^^^^'^'^^l,^,„,^je,u), 6, = ±, ^ e,, = 0. (8) 

i=l 



""The rationale behind this is explained in [23], 
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The two particles form factor of the energy momentum tensor read exphcitly: 

c 27riV^ cosh^ i2L _i2L 

fRSOsi^u) = -^fpiOu) ■ . 1,. \ A e'^'si ® S2 + e ^p^i ® S2), (9) 

smhi(i7r - 6'i2) 

ftsosiOu) = -^/p(gi2)e^(^-+^-) . , (e^.i C5 ."2 + e-i.", s,). (10) 

S smh-(i7r - 6*12) 



and the normahzation for two particles is chosen to be A^® = |iM^Cp, in order to have: 

fnsosi^i + i^' ^1) = 27rM2(ei si (g) S2 + e"^ si (g) 52)- 

One can obtain the form factors of the primaries using the identification s ~ e^^^T^f^'^sai^) ^ 
(thus using the corresponding p- function) then by twisting the Bethe ansatz state like in (jH)) [23] . 



2 Massless perturbation of the coset models 

We consider the massless flows [14] from the UV coset model [15] su{2)k+i su{2)k/su{2)2k+i, 
with central charge 

3A;(A; + l)(2A; + 5) 
''''^ ~ {k + 2){k + 3){2k + 3)' 

to the IR coset su{2)}^ (g su{2)i/ su{2)kj^i. The latter model is the minimal model Mfc_|_2 with 
central charge 

_ k{k + b) 
'''^ " (A; + 2)(A; + 3)- 

The flow is defined in UV by the relevant operator of conformal dimension A = A = 1— 2/(2/c+3); 
it arrives in the IR along the irrelevant operator TT. 

In the massless case, the dispersion relations read {p^-,p^) = ^(e^,e^) for right movers and 
{p^iP^) = ^(e~^', — e~^') for left movers, where M is some mass-scale in the theory, and 6,6' 
the rapidity variables. Zero momentum corresponds to 6 ^ —00 for right movers and 6' — > +00 
for left movers. 

The /S-matrices for the three different scatterings were found in [25] : the RR and LL 5-matrices 
describe the IR CFT and are thus given by the RSOS restriction of the Sine-Gordon S'-matrix. 
The parameter p is related to k hy p = k + 2. The RL scattering is given by^: 



The minimal form-factor in the RL channel satisfies the relation: 

fRL{6 - 6') = fRL{e -e' + 2m)SRL{6 - 
and its explicit expression is given by 

ivr) /■+°° dtl-cosht(l 



t 2 sinh t cosh | 




*For the particular cases p — i and p — +00, the scattering datas were first proposed in [7] and [16] respectively. 
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Its asymptotic behaviour in the infrared is: fniid — ~ /C ea^^ ^' ^'^\ where 

r- I rdt ( 1 

A. = exp — - / — 

2 7o ^ \ sinh t C( 



; cosh I t 

2.1 Form factors of the trace operator © 
2.1.1 Trace operator in the flow from TIM to IM 

The case p = 3 corresponds to the massless flow with spontaneously broken supersymmetry 
between the tricritical Ising model (c = 7/10) and the Ising model (c = 1/2). The form factors 
of have non vanishing matrix elements on an even number of right and left particles. They 
satisfy the residue equation in the RR channel at Oi = 92r + itt {Srr = Sll = ~1)- 

resF2® 2i(^l, • ■ • , Qlr', d'l,--- : 02l) = 



-2ii^2. 



A similar relation holds in the LL channel. The first non vanishing matrix clement has two right 
and two left particles, and reads [8]: 

Ff^,iei,92;e[,e'2) = ^ sinh ^ sinh ^ n /«^(^^ - ^^)- (11) 

i=l,2, 
J=l,2 

The form factors of the trace operator were first computed for a small number of particles in [8] 
in terms of symmetric polynomials. A formula valid for an arbitrary number of particles was 
then presented in [17]. Subsequently in [27] it was noticed that it is possible to simplify the 
formula of [17] such that the form factors could be rewritten as: 

^2r,2l{(^l^ ■ ■ ■ ) ^2r; O'l, Q21) = ^4 



X n si^^v n sinh^n/«^(^^-^i)^2r-,2;(^i, 



2 -L-i- 2 

l<i<j<2r l<i<j<2l i,j 



hr',01, ■ ■ ■ ,0'2l) 



Explicitly: 



Q2r,2«(^l) • • • )^2r; O'l, • • • , 6*2 



21) 



TCS, T'CS', iST^ iST', ieT, ieT'_, 

#T=r--l ^7-/=;_i feeT fegf' keT' keT 

and 

-1 2i . -1 2i 



fRRi0ij)fRRi0ij + i7r) sinh% ' ^^'^ /ll l^^,- ) /ll (^^,- + i^r) sinh 6.. 



r> (12) 



as well as: 
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Let us note that the normahzation of the form factors of the trace operator in (jllj) was chosen 
in order to match the effective Lagrangian [8]: 

where /J and /J^ are the two particles form-factors of the right and left components of the 
energy momentum tensor in the thermal Ising model: 

g (01,92) = 2i7rM2e^^+^^ sinh ^, gie[,e',) = 2i^M^e-''^-''^ sinh ^. (14) 

This normalization for the trace operator ensured a proper numerical check of the c-theorem 
in [8]. 

2.1.2 Generalizations 

The form factors of the trace operator for p arbitrary were first obtained in [17] by Mejean and 
Smirnov in a different representation from the one presented below. The problem of normaliza- 
tion of form factors is not discussed in [17], and the numerical checks are not performed. 
The form factors of G have non vanishing matrix elements on an even number of right and left 
particles. They satisfy the residue equation in the RR channel at 6i = 92r + Itt: 

resF.® .a' a' 



2r,2ilf^l' ■ ■ ■ )f^2r,t7i, • • • ,^21) - 

/ 2r-l 21 
2i F^r-2,2li02, ■■■ , d2r-i;0[, • • • , 621) 1 " Sp^^^{9i - ^2r) SRL{02r - 0'}.) \ S 



R 



1=2 k=l 



p — _ — — _ 

where e = e^psi S2r + e ^^si ^ S2r- A similar relation holds in the LL channel. 

The lowest form factor for the trace operator contains two right and two left particles; it reads: 

-f^^2(^i> ^2; (^1,(^2) 

levrM^ „ T-r „ cosh ^ cosh % (g) 



p2 



YT fpi9l2)fp{0'i2) n fRLiOi - O'j) . ■ N Itai ^ 

^ iti 2 smhi((9i - 6'2 - ivr)smhi(6'i - - ivr) 

n /^^(^^ - e''^^''^~'^''^ flsosiei2)flsos{e[2l (15) 



1=1,2 
J = l,2 



where fRsosi^'^'^) defined in equation (jH 
We make the following ansatz for the solution of the form factors equations with the first 
recursion step given by (fT5|l : 

^2®, 2/(^1 > ■■■,(^2r]0[,...,92i)= jj'2^^4: H /p(%) U /p(^ij)Il//JL(^i - ^j) 

r I 
JCg j^^i JCg/ 



9 



The p-functions P^5os(^''") defined through equation and the analogue of equation 0. 
The Bethe ansatz state ^^{9,u) was introduced in (jH)). 

At 01 = + ivr, each of the r integration contours gets pinched at 92r, ^2r i ivr. For this reason 
we introduced the function M.2r,2i{Q-, S', u, v) that should satisfy the properties at 6i = 6*2^ + ivr: 

M2r,2l{0i, . . . ,02r;O';Ui, . . . ,Ur]v) = 

21 



M2r~2,2l (6*2 , • • • , 02r-i; 6';Ui,... , Ur-l] v) 

• Ur = 92r ± ivr: 



k=l 



M2r,2l{0l ■ ■ ■ 02r; d'] Ui . . . Ur] v) = 

M2r-2,2l{G2, ■ ■ ■ , ^2r-i; Ul, - ■ ■ , Ur-l',v) 



21 



k=l 

Similar relations hold in the LL channel. 

We introduce the sets S = {1, . . . , 2r) and 5' = (1, . . . , 2/), as well as the subsets T C S ,T' C S' 
and T = S\T, T' = S'\T' . These subsets have the number of elements: #T = r — 1, #T = r + 1 

T = {ii < ^2 < • • • < V-i}, f = {ki < k2 < ■ ■ ■ < kr+i}, 

and T',T' are defined similarly. We propose the following function: 

A^r+l-2 n^.^efCOS^ UkAeT'COS^ 

M2r2i{e,e',u,v)= „ , y y Ae^Ee..^<^ — e^S^'L 

#T=r-l #T' = ;-1 kef 

n ^T, cos^^n cos^^ 
X n "'"'"%oc"-"-n'"^'""' CO.---- n m - n ^(^^ - ^^') • 

lll<m<n<r '""^ 2i I ll<m<n<l '^^^ 2i isT, ,gt'. 

In particular, for two right movers and two left movers: A^2,2(^i) ^2! S'l, ^2! ^) ~ e~^^~^^^^i~''^2. 
A similar function was first obtained in [17] in a different representation^, and was used in a 
very different context in [26]. We present it here in what we believe to be a simpler expression. 
It should also be compared with a very resemblant function obtained in [12] in the massive case. 
Finally, the following recursion relations hold [20]: 

Nfr = Ng-2{fr\^)?l^^r, = Ni_,{f^^-{0)f/47rl. 

p, 2ih'PC'^ 

We set the normalization of the p-function N2 = — ^5—^ in order to have in the infrared the 
same relation as in the flow from TIM to IM^'^: 

^2 ^ ~^^-^2^-^2^' 



^The integral representation for the form factors proposed in [17] is different from ours, as it contains two 
integration variables less (one less for each right/left channel). We did not manage to prove that the representation 
for our form factors coincides with the one of [17]. 

^"We are grateful to G. Delfino for explaining to us how to normalize properly the form factors of the trace 
operator. 
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where f2 and f2 are the two particles form-factors of the right and left components of the 
energy momentum tensor in the Mk-\-2 minimal model perturbed by ^1^3 -see (|in|) -. themselves 
normalized such that: 

In [17], it is shown that for an arbitrary number of particles the form factors of the trace operator 
reproduce the direction of the flow: 2« ~^ ~^;^W^ f2rf2i- verified with Mathematica that 
this is true on our representation, but this check was done for a small number of particles only. 

3 Numerical results 

The knowledge of the form factors of the trace of the stress energy tensor allows to estimate the 
variation of the central charge along the flow by means of the "c-theorem" sum rule [5,6]: 

o f-CO 

Ac = Cuv-CiR = - / drr3(e(r)e(0)). (17) 

Since in the massless case any correlation function can be represented by its spectral expansion 
(O(x)O(O)) = (18) 



00 



r,l=0 ^ 

the computation of Ac turns out to be a non trivial check for the form factors F^2i (^t least 
for the first few of them) . 

In the past years, the calculation of Ac using form factors has been done in various massive 
theories, providing accurate results even within the two-particle approximation. However in [4], 
and more recently in [12] within the construction of form factors for the (massive) SS model 
and its RSOS restrictions, significantly large discrepancies have been observed comparing the 
computation of the central charge by means of the c-theorem and the corresponding exact results. 
In [12], within the two-particle approximation and for a given subset of the parameters of the 
model, the deviations were of about 20-25%. On the other hand, the massless flow from the 
Tricritical Ising to the Ising model provided up to now the unique massless case known where 
numerical checks for the c-theorem were performed: the numerical results obtained in [8] by 
Delfino, Mussardo and Simonetti are unexpectedly accurate, as the leading contribution (four 
particles) was enough to obtain the 98% of the exact result. This is a very remarkable situation, 
as there is absolutely no reason to expect in the massless case that the leading contribution 
gives any good approximation at all to the correlation function: at a given energy, one cannot 
say what number of particle processes contribute. 

In the present case, we deal with a one parameter family of massless flows, whose variation 
of the central charge is given by 

(we recall that k = 1 corresponds to the flow from TIM to IM). Moreover, it is of interest to 
compare the accuracy of the numerical results for the variation of the central charge in the 
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Table 1: Massless flow - four particles approximation (the number in bracket in the second 
column is the statistical error on the previous digit due to the Montecarlo integration). 



massless case against the accuracy obtained in the associated massive coset model [12], where 
the UV CFT is the same as the one considered in the previous massless flows, and whose central 
charge is given by 

_ 3fc(fc + l)(2A: + 5) 
" ''^^ " {k + 2){k + 3)i2k + 3)- 

Let us start with the massless flow: in the following we will consider the correlation function 
truncated to 4 particles, with the use of the form factor (|15|) : higher form factors are very 
difficult to compute numerically because the scattering is non-diagonal for k > 1 (the 6-particles 
contribution for the diagonal case k = 1, corresponding to the TIM — IM flow, was considered 
in [8], starting with the formula (|11() for four particles). 

Within such an approximation we obtained the numerical estimates collected in Tabled where 
the comparison with the expected results (|19|) is also given. In order to explain the observed 
deviations from the exact values of Ac^""""*, it is worth remembering that the conformal dimension 
of the perturbing operator (in the UV) is given by 

9 

A = A = 1 



2^ + 3 



which becomes marginally relevant in the limit /c — > oo. Hence, in such a limit, the UV behaviour 
of the correlator in the sum rule tends to (also logarithmic corrections may appear), making 
the convergence of the integral (|17)) weaker as k is increased. 

As a consequence, at large but finite k the suppression in the UV region (i.e. where the 
form factors approximation is worse) becomes very weak giving rise to the observed deviations. 
In principle, such an argument can also explain the reason why the four-particle approximation 
shows large discrepancies at relatively small values of k. 

In order to illustrate this point let us consider the following example: when k = 4 we have 
A = 0.818 . . . and a corresponding deviation of the 12% which can be compared with the case 
k = 1, where the deviation is of about 2% and one has A = 0.6. Hence it is likely to conjecture 



12 



that the previous mechanism is a good candidate to be the responsible for the poor quahty of 
the numerical estimates of Ac^. Such an interpretation also implies that it is of a moderate 
interest to work out the contributions due to 6-particles form factors: we expect that they will 
add a quantitative correction to our results without changing the global picture outlined above. 

A similar behaviour has been observed [12] in the massive flow originating from the coset 
model^-^ su(2)fc+i (g) su(2)fc/sii(2)2fe+i perturbed by the same operator as before with conformal 
dimension A (the sign of the perturbation is different). In the massive case, the lowest form 
factor needed for numerical tests contains 2 particles. As one can see in Table |21 the numerical 
determination of the central charge given by ()2U() becomes worse as k is increased (we observe 
that in this case the discrepancy is not as large as in the massless flow^^). Since the perturbing 
operator is the same, we expect that the reason of such a decreasing of the precision is due to 
the same mechanism as before. 

As a side remark it is interesting to notice that the limit k ^ oo seems to resemble the 
case of the SU{2) Thirring model (which is nothing but the Sine-Gordon model at the marginal 
point /?2 = Svr): even if the perturbation is marginally relevant, the two-particle approximation 
to the sum rule gives a finite result. The very important difference is that in the latter case the 
numerical computation [12, 28] gives c*-^^ = 0.987 which is unexpectedly near to the exact one 
goxact _ ^ (instead of being very far as in the case of coset models, where we have 43% and 21% 
for the massless and massive models respectively). The previous examples shows that one has 
to be very cautious in interpretating the results which come from the form factor approximation 
to the c-theorem sum rule in the case of marginal perturbations (we refer the reader to [12] for 
additional numerical examples in the S'C/(2)fc-Thirring model)^^. 

We find it interesting to present the results for the case k = —1, which corresponds this time 
to a non-unitary flow from a UV fixed point with Cuv = to a IR CFT with Cir = —2; the 
result in Table Q shows that the accuracy of the estimate for the central charge is good, and 
the fact that the approximation to the exact value comes from above is in agreement with the 
non-unitarity of the flow. 

Finally, we used the 4-particle form factors to compute in the massless case the correlation 
function (0(x)0(O)). In the flgures 1-4 we compared such an approximation of the correlator 
with its power-law behaviour at both the UV and IR flxed points for different values of A; (A; = 1, 
2, 3, 10; the case k = 1 was already considered in [8]). One may observe that the relation (fTC]) 
plugged into the integral (|18|) followed by a simple change of variables insures the IR behaviour 
of the correlation function taking into account the four particles contribution only, 



^^This coset model is obtained by RSOS restriction of the SS model [13]. 

^^Such a difference is due to the presence of the exponential factor exp(— mi? cosh 9) inside the spectral expansion 
for correlators in massive theories, see [8]. 

^^More generally, even within the (massive) minimal models Mk+2 one observes a loss of precision when in- 
creasing the parameter k [12]; however it is nothing comparable with the phenomenon observed here. 




where 
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Table 2: Massive flow - two particles approximation 



As for the UV behaviour, it is given by the following conformal OPE: 



Ok (Mr) 2fc+3 . 



(e(x)e(o)) ~ M 

r->o \ 2k + 3 

where the constants can be extracted from [29] (Eq. (2.2) with h = n = 2 and I = k + 1) 
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Looking at the diagrams 1-4 one can clearly notice the difference between the expected 
UV power-law behaviour and the slope given by the truncated form factor expansion. As 
a consequence, the combined effect of such a difference in the UV behaviour and the above 
mentioned weak suppression in the integral (|17|) can be considered as the cause of the discrepancy 
between c-theorem and exact results. 

Remarkably enough, it appears on the charts that at the condition of not going too far in 
the UV, the four particles truncation seems to give a reasonably reliable approximation to the 
correlation function all the way from the intermediate region to the (exact) IR region. 



Concluding remarks 

In this work we constructed form factors of the trace operator in the massless flow between 
the coset model ,3^(2)^+1 ^ su{2)k/su{2)2k+i and the minimal model. We mimicked the 

construction of form factors of the trace operator in the massless flow TIM — s- IM obtained in [8]: 
in the section 2.1.2, we looked for a solution of the residue equation obtained by replacing the 
RR and LL ^-matrices S^^^^ = — 1 by the 5-matrix S^^2^, that should match the desired 
properties in the IR. The form factor of the trace of the stress-energy tensor Q with four particles 
was used to compute both the (truncated) correlation function (B(x)0(O)) and the variation of 
the central charge (by means of the "c-theorem" sum rule). Since for the latter we observed a 
large discrepancy with respect to the exact results, we carefully analyzed the problem, giving 
a plausible explanation for the phenomenon. Nevertheless, it appears on the diagrams 1-4 that 



14 



the truncation to 4 particles for (0(x)0(O)) still gives not too bad an approximation to the full 

correlation function, at the condition that one does not go too far in the UV. 

In our next article [30], we will generalize the construction of form factors of the magnetization 

operator and the energy operator in the massless flow TIM IM [8] to the whole family of 

flows. 
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Figure 1: Flow k = 1, logarithmic plot of the correlator G{r) = {&{x)@{0)) (black line) together 
with both the IR (grey line) and the UV (dark-grey line) behaviours. 
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Figure 2: Flow k = 2, logarithmic plot of the correlator G{r) = (G(x)0(O)) (black line) together 
with both the IR (grey line) and the UV (dark-grey line) behaviours. 
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Figure 3: Flow k = 3, logarithmic plot of the correlator G{r) = {&{x)@{0)) (black line) together 
with both the IR (grey line) and the UV (dark-grey line) behaviours. 
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Figure 4: Flow k = 10, logarithmic plot of the correlator G{r) = {Q(x)Q{0)) (black line) together 
with both the IR (grey line) and the UV (dark-grey line) behaviours. 
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